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Abstract 

In this paper we treat the Jaynes-Cummings model with dissipation and give an ap- 
proximate solution to the master equation for the density operator under the general 
setting by making use of the Zassenhaus expansion. 

In this paper we treat the Jaynes-Cummings model (pQ) with dissipation (or the quantum 
damped Jaynes-Cummings model in our terminology) and study the structure of general 
solution from a mathematical point of view in order to apply it to Quantum Computation 
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or Quantum Control, which are our final target. As a general introduction to these topics 
see for example |2] and |3]. 

We believe that the model will become a good starting point to study more sophisticated 
models (with dissipation) in a near future, see for example Concluding Remarks in the 
paper. 

Let us start with the following phenomenological master equation for the density operator 
of the atom-cavity system in |1] : 

= —i[Hjc, p] + p |apa^ — -(a^ap + pa'''a) j> + u i^a^pa — -{aa)p + paa^) | (1) 

where Hjc is the well-known Jaynes-Cummings Hamiltonian (see |T]) given by 

Hjc = yds ® 1 + u;ol2 ® a"^a + ^^ (cT+ ® a + cT_ ® a^) (2) 
f + LOoN na 



with 



oi\ 

ooy V^O/ 

, and a and are the annihilation and creation operators of an electro-magnetic field mode 
in a cavity, = a'^a is the number operator, and p and z/ (p > z/ > 0) are some constants 
depending on it (for example, a damping rate of the cavity mode). 

Note that the density operator p is in M(2; C) ® M{J^) = M(2; M{J^)), namely 

p^f"" ''°'),M(2;M(^)) (3) 

\ PlO Pll / 

where M(J^) is the set of all operators on the Fock space J-" defined by 



J- = Vectc{|0),|l),|2),|3),---} 

oo oo 



. n=0 n=0 

and 1 is the identity operator 



We would like to solve ([T]) explicitly. However, it is very hard to solve at the present 
time, so we must satisfy by giving some approximate solution to it. 

First of all we write down the equation ([1]) in a component-wise manner. For that we set 



H 



Qa 



JC 



fiat -^ + uoN 



A B 
C D 



for simplicity. 

Then it is easy to see 

Poo = -^(^Poo + BpiQ - pqqA - poiC) + 

P japoofl^ - ^(a"^apoo + Pooa"^a)| + jaVooa - ^(aaVoo + Pooaa^)| , 
: -iiApoi + Bpu - pooB - poiD) + 



Poi 



Pio 



Pii 



1 



/i < apQio) — -{a^apoi + pQia)a) \ + { Poia — -(aa^poi + poicto-^) 



1 



-^(Cpoo + Dpio - piqA - pnC) + 

p i^apioa^ — -{a^apiQ + piQa^a)^ + z/ jatp^ga — -(aaVio + Pio^a^) 
~i{Cpoi + Dpn - piqB - puD) + 

p lapiio) — -{a)apii + pno^a) 1 + iatpua — -{aa) pu + puaa)) 



(4) 



^ (5) 

where pij = {d/dt)pij as usual. 

Here let us remind some technique used in |2] and which is very useful in some case. 
For a matrix X = (xij) G M{J^) 

( \ 

Xoo a^oi 2;o2 ■ ■ ■ 



X 



Xio X\\ X\2 ■ ■ ■ 
X20 2^21 2:22 ■ ■ ■ 



we correspond to the vector X G J^dimc^^ g^g 



^ — (a^ij) — ^ — (xoo, a^oi, a;o2, ■■■ j a^io, a;ii, X12, 2:20, a;2i, 2:22, ■■■ j ^ (6) 
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where T means the transpose. Then the following formula 



EXF = ® F^)X 



(7) 



holds for E,F,X E M{J^). 
This and equations ([5]) give 



/i <j a (g) (a"^)^ - ^{a^a (g) 1 + 1 (g) a^a) \ + u la^ ^ - ^(aa^ (g 1 + 1 (g aa^) 



PoO) 



poi = B^poo+ D^) Poi + B ®lpu} + 

/i I a (g (a"^)^ - ^(a"^a (g 1 + 1 (g a"^a) \ + u ^a^ - ^(00"^ (g 1 + 1 (g aa"^) 



Poi, 



^10 = -i{C®lpoo+ {D(^1-1(^A'^) Pio-l®C^Pu} + 

/i la (g (a"^)'^ - l-{a^a (g 1 + 1 (g a"^a) I + z/ | a"^ (g - ^(aa"^ (g 1 + 1 (g aa"^) 



ySii = -i{C ^lpoi-1® B^pio+ {D ^1 -1^ D^) pii} + 

/i I a (g (a^")^ - lr{a^a (g 1 + 1 (g a"^a) I + z/ | a"^ (g - ^(aa"^ (g 1 + 1 (g aa^) 



Pii 



(8) 



because 1 and = a'^a are diagonal (1"^ = 1, N'^ = N). 
From ([3]) we set 

Poo 
Poi 

PlO 

ypii / 

then we obtain the following "canonical" form 



Poo Poi 



Pio Pii 



(9) 
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+ 



A^l-l^A^ -lOC^ B^l 

AOI-IOD^ 5(8)1 



(10) 



with 

L = ^ |a (g) {a'')^ - ^(a"^a (g) 1 + 1 (g) a'^a)! + u i^a^ ^ - ^{aa^ 1 + 1® aa') 

or more explicitly 

d „ 



f oja{N ^ 1 ~ 1(g) N) (g) (flt)^ f2a(g)l N 

-ni(S)a^ ojQ + ojo{N (g) 1 - 1 (g) N) fia (g) 1 

na'^(g)l -cjo + wo(^® 1 - 1 ® ^) -ill (X) (a^')^ 

\ fiat 1 -fil(g)a^ a;o(iV ® 1 - 1 (g) iV) / 

/ L \ 



+ 



from 



L 



(11) 



L J 



/I = y + wqA^, B = na, C = Qa\ D = -y + coqN. 



Here a simplified notation uq in place of WqI (S) 1 has been used. 
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Now we set 

( loq{N® 1-1® N) (g) (flt)'^ na®l \ 

-Vll®a^ ijjQ + ljJq{N ® 1 - 1 ® N) Vla®l 

Vta^ ®1 -loq + ujfi{N ® 1 - 1 ® N) -m®{a^Y 

®1 -m®aF llIq{N ® 1 - 1 ® N) J 



X 



Y 



\ 

( L \ 

L 

L 

V L J 



(12) 



(13) 



for simplicity. 

Next, in order to rewrite X and Y in terms of Lie algebraic notations used in [5] we set 



aJ®a\ K_ = a®{a^f, fsTs = -(iV ® 1 + 1 ® iV + 1 ® 1), 



(14) 



Then it is easy to see 



(ir+)t = ir_, {Ksy = K,, {Koy = Ko, 
[K,,K+]=K+, [K,,K^] = -K^, [K+,K. 
[Ko,K+] = [Ko,K^] = [Ko,Ks] = 



(15) 



and 



[/sTo, a ® 1] + a ® 1 = 0, [Kq, 1 ® a^] + 1 O = 0, 

[Ko, ® 1] - ® 1 = 0, [i^o, 1 ® (a^)^] - 1 ® (a^)^ = 0. (16) 

Namely, {K^, K^} are generators of the Lie algebra S'u(l,l), see for example [7] as a 
general introduction. 

Under the above notations X and Y in (1121) and (IT^ can be rewritten as 
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X 



and 



/ 



Y 



UqKq -Q1 (g) (a"f)^ Qa^l 
Vta) ®1 -Vtl (g) ojqKq 



\ 



\ 



(17) 



/i — z/ 



+ {uK^ + — (/i + z/)K3) I4 (for simplicity) 



:i8) 



where in the process the relation aa^ = a'' a + 1 = + 1 has been used. 
As a result we obtain the final form 

d „ 



dt 



P 



-ix + yYp 



(19) 



with X in 0171) and Y in fllSp . This is our main result and we believe that the equation is 
clear-cut enough (compare it with ([1])). 
Since the general solution is given by 



p(t) = e*(-'^+^)p(0) 



(20) 



in a formal way we must calculate the term e**^ ^x+y) ^ which is in general very hard. For 
that the following Zassenhaus formula is convenient. 



Zassenhaus Formula We have an expansion 



-V{2[[A,i?],i?]+[[A,i?],A]}^V[^.^]e*^e*^ 



(21) 



The formula is a bit different from that of [8]. 



In this paper we use the approximation 



with A = —iX and -B = F, so that p(t) is approximated as 
Let us calculate each term explicitly : 



(22) 



[I] First, we calculate e However, the calculation is more or less well-known. We 
decompose X into two parts 

\ 



—VLI ® ojq + ujqKq 























X,. 








na^i 

Qa^l 



\ 



Vta^ ®1 
fiat (g) 1 









Then it is not difficult to check 



owing to the relations ( jTBjl . If we set 



[Xi,X2] = 



s 



10 
10 
10 
1 



then 



sxoS'^ = sxoS = n 







a® 1 
at (g) 1 

at ® 1 





a (g) 1 



(23) 



(24) 



(25) 



(26) 
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Moreover, if we decompose like 

-fll <^ ojq + uqKq 
then 



UqKq 

Uo + UqKo 



+ 



-Vtl 



-Ql ® {aY 
ni (g) 



UqKq 

Wo + wo-fCo 

owing to the relations (1161) . Therefore we have a decomposition 



n 



\ 



uoKo 

LOo + UoKo 

-ujQ + loqKo 

iOoKo 

l(8)(at)^ 

l®a^ 

l^ia^f 

l®a^ 



As a result we have a decomposition consisting of three commutative operators 



X 



uoKo 

Uo + LOoKo 

-ooo + u}qKo 

ujqKq 

l®{a)Y 

l®a^ 

l®{a)Y 

l®a^ 



\ 



+ 



n s 



a®l 
a+Ol 

a®l 

a) ®1 
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Now, by making use of the well-known formula 

A \ \ / cos{aVAA^) -i^j^sm{aVAA'<)A 

-1-7= sm{a^/A^A)A^ 



exp —la 



A^ 



cos(q;\/AU) 



where a is some parameter, we can calculate e easily. The result is 

^(11) (12) (13) (14)^ 



-.-itX 



(21) (22) (23) (24) 
(31) (32) (33) (34) 
(41) (42) (43) (44) ) 



where for a = VLt 



(12) 



(13) 



^-ituioKo cos(q; V'T^oot) cos(a\/aat 1) 

^-itu^oN cos(q;V]vTT) ® 6**'^°^ cos{a^^WTl), 

■ ^-itu,oKo_ ^ sin(a\/l (g)aat)(l ® {a^f) cos(a\/aat (g) 1) 



= i 6-'*'^°^ cos(q;VA^ + 1) ® e'*'"° 



AT 



1 



VN + T 



sm{aVN + l){a'f, 



-i e~^*^°^°cos(QiVl (g) ggt) ^ sin{a^/ aa^ ® l){a ® 1) 

V aat (8) 1 



= —le 



-itwoN 



sm{aVN + l)a ® e^*^"^ cos(aVA^+ 1), 



(14) = e-^*'^°'^°^^=sin(Q;Vl (g) aat)(l ® (a^)^)^^= sin(Q;A/aat (g) l)(c 
- i : sm{ay/NTl)a (g) e^*'^"^ , ^ -. sm{ay/NTl){aY; 



and 
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(21) = ie-^^Oe-'^°^° ^ sin(aVl®ata)(l(8)a^)cos(aVaat0 1) 
= ie-^*'"° 6-'"^°^ cos(q;\/]vTT) (8) e'*'^°^^ sin(Q;y]V)a^, 



(22) = g-iia;og-itu;oXo cos(Qi Vl ® a^a) cos{a^/aa^ 1) 
^ g-«ta,o g-itcoiv cos(q;\/]vTT) ® 6'"^°^ cos(q;\/]V) , 

(23) = e-^^^Qe"^*"^"-^" sin(a\/l ata)(l (g) g^) — sm{a^/ aa^ ^ 1) {a ^ 1) 

= e-^*'"° g-ita;oiv^^^ sin(a\/iVTT)a ® e**^°^^ sin(a^)a^, 

1 



(24) = _i^-itioo^-itu,oKo cos^q;^! (g) atg)^ sm{a^/aa^ l)(a (g) 1) 

vaat (g) 1 

= -ie-^*'"^ g-ita;oiv^^^ sin(Q;V']vTT)a (g 6^*^^°^ cos(q;V]V); 

and 

(31) = -ie^*'^°e~^*"°-^°cos(QiVl0aat) , — sin(Q\/ata g) l)(at g) 1) 

Vata (g 1 

= -ie'*"° e-^*'"°^^ sin(a ViV)a^ ® e**""^ cos(a\/iVTT), 



AT 

(32) = g»toog-itooj^o ^ sin(Q;A/l(gaat)(l(g(at)^)— sm{a^/a^a l){a^ <^ 1) 

V 1 (g aat Va'I'a (g 1 

= 6^*"^° e-i*a;oiv 1 sin(Q;V]V)at (g e'*^°^ , : sm(ay/NTl){aY , 



(33) = e'*'^°e-^*'^°^° cos(a\/r¥aat) cos(a\/at^^) 
= e^*'^" 6"^*'"°^ cos(q; VtV) (g 6^*^^°^ cos(q; V^VTT) , 

(34) = jgitoog-iMo 1 sm(aVl0aat)(l0(at)^)cos(aVata(gl) 

V 1 ® oa^ 

= ie^*'^° 6-^*'"°^ cos(q;V]V) (g g'*'^"^— ^= sin(Q;v/]vTT)(a^)^; 

ViVH- 1 

and 
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(41) = e~^^'^°^°^^=sin(aVl^ata)(l0a^) ^ = sin(a Vata l)(at ^ 1) 
V 1 (g) a^a Va'I'a 1 



(42) 



(43) 



-UujqN 



sin(Q;v N)a: ® e 



sin(av N)a 



e 



-*ta;oiv 1 sin(ayiV)at ® e**"«^ cos(aViV), 



AT 



■^-itujoKo sin(a\/l ata)(l ® a^) cos(aVata ® 1) 

V 1 ® a^a 



sin(av N)a 



(44) = e-**'^«^o cos(a\/l ® ata) cos(a v^ata O 1) 
= e-^*"^"^ cos(ayiV) ® e^*"'"^ cos(ayiV). 



[II] Second, we must calculate e*^. To be glad this task has been done, see [5] and [6]. 
Namely, from ffTSjl we have 



and 



\ 



\ 



4L 



L = ^^—^ + vK+ + ^iK_ - (/i + z/)K3 (31) 



^ M^t^GWA-+ g-2 \og{F{t))Ks ^E(t)K- 



(32) 



with 



E(t) 

m 

Git) 



-^sinh(^t) 



cosh(i^t)+g^sinh(^t) 
cosh ( : — 1 \ + — sinh | : — t 



2 J fi — u 
^sinh(^i) 



cosh Sinh (^t) 

12 



(33) 



This is a kind of disentangling formula, see for example [7] as a general introduction. 



[III] Third, we must calculate e"*^'"^'^^ In this stage some interaction appears. Let us 
calculate [X, Y] exactly. 
Some calculation gives 



[1 (g) a^, vK+ + - (/i + y)K:i] = -fia (g) 1 + 



1(g) 



1 ® (a^) 



[a (g) 1, vK^ + - (/i + z/)ir3] = ul®a^- 



a (g) 1, 



(34) 



and for simplicity we set 



A = -Z/flT (g) 1 + ® (flT)^, 

B = fxaig)!-- l(g)a , 

C = z/1 (g) a - ^ a (g) 1, 

D = -;,l®(atf + ^at®l. 

Note that i? 7^ —A"'' and D ^ —C^ because of > z/. It is easy to see 



(35) 



[A,C] = [A,D]=0, [B,C] = [B,D] = Q, 



[A, B] = [C, D] 



(36) 



From the decomposition fl28|) we have a decomposition consisting of two commutative 
operators 



[X, Y] = n 



A ^ 

5 

A 





c ^ 

L> 

000c 

L> 



with ^ in (ESI). 



(37) 
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Now, by making use of the formula 

for some parameter /3, we can calculate the term e~*~['^'^l easily. The result is 

^(11) (12) (13) (14) ^ 



where ior /3 — 



(21) (22) (23) (24) 
(31) (32) (33) (34) 
\^ (41) (42) (43) (44) J 



(11) = cos{(3^/AB) cosipVCD), 

(12) = -i--^sin(/3\/AB)^cos(;5A/aD), 

(13) = -icos(^VAB)-^i=sin(^v^)C, 

(14) = --^sin(/3VAB)^-^ sin(/3 VaD)C ; 

(21) = -i-^sm{P^/BA)B cos{/3VCD), 

(22) = cos(/3VbA) cos(P VCD), 



(23) = --^sm{P^/BA)B-^smif3VCD)C, 

(24) = -icos{pVBA)-^= sm{^ VCD)C ; 



(31) = -icos(^VAB)^^sin(^v^)L>, 

V DC 

(32) = --^sm{l3VAB)A-^sm{l3VDC)D, 

(33) = cos(^VAB) cos(P^/DC), 

(34) = -i-^sin(/3VAB)^cos(/3/DC) ; 
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(41) = 7^sin(/3/BA)5— ^sin(/3/DC')/^, 

Vba Vdc 

(42) = -icos{f3VBA)^^sm{f3VDC)D, 

V DC 

(43) = -t-^sm{P VbA)B cos{PVdC), 

(44) = cos(/3/bI) cos(/3/DC) 



where 

y4i? = -^va)a (g) 1 + ^ — z/a^ (g) + /i^—^ — a ® {a^y - ( ^ — j 1 (g) aa^ 
= -iiuaa^ (g) 1 + /i^-^ — a (g) (a^) + ^-^ — ua^ (g) a - ( ^—^ — 1 1 (g) a^a, 

2 



CD = (g) a^a + ^— ^ — z/aT (g) + ^ — a (g) (a^)-* - ( ^— ^ — j aa^ (g) 1, 

^ + « + , M + z^ + T fu + v\'^ + 

L'C = -/iz/1 (g) aa} + a (g) (a^)^ + ^— ^z/a^ (g) a - ( j a^a (g) 1 

These are comphcated enough. 

Anyway, we have completed the task although it is very complicated. 

In last, we shall restore the result to original form. If we set 

^(f)=e*^e-*^p(0) = e*^p\(t), 
then pi{t) is from fl30|) given by 



(11) (12) 
(21) (22) 



where 
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(11) = g-it-'oAf cos(a ViV + 1) Poo e**"^"^ cos{aVN + 1) 

+ie-'*"«^cos(ayiVTl) Poi ate'*"»^^^= sinraViVTT) 

ViV + 1 



and p(t) is by 



-UujqN 



sin(aVA^ + l)a pio 6^*^°^ cos(aVAr + 1) 



+6-^*'^°^ —^ sin(a ViVTT)a Pn a^e'*'^"^ — ^ sin(av/iVTT), 

(12) = ie"**"^" e-^*'^o^ cos(a VAT + 1) poo ae^*"""^^ sm(a\/iV) 

V N 



+e-**"o e-'*"°^cos(ayiVTT) Poi e'*"«^ cos(ayiV) 

+e-'*"» e-^*"°^ — ^ sm(aViVTT)a pio ae'*"«^^ sin(ayiV) 



, ^ sin(aViVTT)a pn e**'^"^ cos(ayiV), 



(21) = -ze**'^« e-**'^«^^sin(ayiV)a"f Poo e'*'^°^cos(aViVTT) 



Vn ^/NTT 



+e''^° e-'*"o^cos(ayiV) pio e'*"°^ cos(ayiVTT) 

+ie^*'^o g-ita;oiv cos(ayiV) Pn a'fe^*'^°^^^= sm(ayiVTT) 

(22) = e-'*'"o^^sin(ayiV)at Poo ae^*"'o^^sm(ayiV) 



.^e-^*-o^^sin(ayiV)at p^^ 6^*"°^ cos(ayiV) 
v A^ 



+ie-'*""^ cos(ayiV) pio ae'*"o^-^ sin(ayiV) 



^g-*ta;oJv cos(aVAr) Pn e**"^°'^ cos(aVAr) 



^ n=0 ■ lm=0 ■ J 

exp({-log(F(t))}A^)]a" (42) 

in terms of E{t), F{t), G{t) in fl5^ and A^ = a'''a. See [3] and [B]. This form is relatively clear 
because of no interaction. 
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Our (approximate) solution is 

from fBU]) . We would like to express p{t) like f l42l) . From fl5^ we can expand each term in 
terms of the Taylor expansion of cos(X) and sin(X). However, the form is very ugly in the 
Dirac's sense, so we leave such a task to readers. 



Concluding Remarks In this paper we treated the master equation for the density op- 
erator based on the Jaynes-Cummings Hamiltonian with dissipation and constructed the 
approximate solution up to 0{t^) under the general setting. Our construction is quite gen- 
eral because p(0) is any initial state. 

However, it may be still inconvenient to apply to realistic problems coming from the 
atom-cavity system. Further work will be required. 

Moreover, it may be possible to treat a generalized master equation for the density 
operator given by 



—p = —i[Hjc,p]+l^lapa^ — -{a^ap + pa^a)\+i'<a''pa — -{aa''p + paa^))- + 



n\apa- -{a^p + pa^) j + k ^a) pa) - -{{a^fp + p{,a)f) | (43) 

with the positivity condition p,v > However, we don't treat this general case in the 

paper, see [9], [10] and [TT] . 
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